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Abstract 

We study the decay width, forward-backward asymmetry and the longitudinal lepton 
polarization for the exclusive decay B — ► K*l + l~ in the two Higgs doublet model with 
three level flavor changing neutral currents (model III) and analyse the dependencies of 
these quantities on the the selected parameters, £ U,D , of model III including the next to 
leading QCD corrections. It is found that to look for charged Higgs effects, the measure- 
ments of the branching ratio, forward-backward asymmetry and the longitudinal lepton 
polarization for the decay B — > K*l + l~ are promising. 
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1 Introduction 



Rare B meson decays induced by flavor changing neutral current (FCNC) b — > s transition 
represent an important channel to test the Standard model (SM) at loop level. They provide 
a comprehensive information in the determination of the fundamental parameters, such as 
Cabbibo-Kobayashi-Maskawa (CKM) matrix elements, leptonic decay constants, etc. and open 
a window to investigate the physics beyond the SM, such as two Higgs Doublet model (2HDM), 
Minimal Supersymmetric extension of the SM (MSSM) U, etc. Measurement of the Branching 
ratios (Br) of the inclusive B — > X s 7 [0 and the exclusive B — > K* / y |3[] decays stimulated the 
study of the rare B meson decays in a new force. 

Currently, the main interest is focused on the rare decays where the SM predicts large Br, 
which is measurable in the near future. B —>■ K*l + l~ decay is one of the candidate of these 



decays and it is described by b — > sl + P transition at the quark level. In the literature [Q- JT7 
this transition has been investigated extensively in the SM, 2HDM and MSSM. The forward- 
backward asymmetry of dileptons is another quantity which provides information on the short 
distance contributions. Further, the longitudinal lepton polarization includes different combi- 
nations of the Wilson coefficients C% ,Cg and Ciq and it can shed light on the investigation 
of their signs in the SM and new physics beyond. 

The theoretical analysis of exclusive decays is more complicated due to the hadronic form 
factors, which brings sustantial uncertainity in the calculations. However, it is well known that 
the experimental investigation of inclusive decays is more difficult compared those of exclusive 
ones. Therefore we will consider the exclusive B — > K*l + l~ in the present work. 

The calculation of the physical observables in the hadronic level needs non-perturbative 
methods to determine the matrix elements of the quark level effective Hamiltonian between 
the hadronic states. This problem has been studied in the framework of different approaches 
such as relativistic quark model by the light-front formalism ||17|| , chiral theory fl8| , three point 
QCD sum rules method fl9 |, effective heavy quark theory [^(J and light cone QCD sum rules 
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In this work, we present the next to leading order (NLO) QCD corrected effective Hamil- 
tonian in the 2HDM with flavor changing neutral currents (model III) for the inclusive decay 
b — > sl + l~ and calculate the differential Br of the exclusive B — > K*l + P decay. Further, we 
study the forward-backward asymmetry (Aps) and the longitudinal lepton polarization (Pl)- 

The paper is organized as follows: In Section 2, we give the NLO QCD corrected Hamiltonian 
and corresponding matrix element for the inclusive b — > 57 decay and calculate the matrix 
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element. Section 3 is devoted to the calculation of the differential Br, A FB and Pl of the 
exclusive B — > K*l + l~ decay . In Section 4 we analyse the dependencies of the Br, Afb and 
Pl on the the Yukawa couplings £^ and estimate the BriB — > K*l + l~) for the model III. 

2 QCD corrected short-distance contributions in the model 
III for the decay b — > sl + l~ with additional long-distance 
effects 

In this section we present the NLO QCD corrections to the inclusive b — > sl + l~ (I = e,fi) 
decay amplitude in the 2HDM with tree level neutral currents (model III). Before presenting 
the details of calculations, we would like to give the essential points of this model. 
The Yukawa interaction in this general case is 



C Y = vYjQiL^iUjR + r/FjQiKpiDjR + gjQiLfcUjR + ZgQiLfaDjR + h.c. 



(1) 



where L and R denote chiral projections L(R) = 1/2(1 =p 75), 0j for i = 1, 2, are the two scalar 



doublets, r]fj D and are the matrices of the Yukawa couplings. The Flavor Changing (FC) 
part of the interaction can be written as 



-U,D 



£>Y,FC = {.ijQiL&UjR + ZtjQiL&DjR + h.C. 



D, 



(2) 



with the choice of 0i and 02 



V2 





v + H° 



+ 



( V2x 



Here the vacuum expectation values are, 



< 0i >= 



1 



(° 



V2 V v , 

and the couplings £ ' for the FC charged interactions are 

£ch = ^.neutral VcKM , 



< 02 >= , 



neutral i 



(3) 



(4) 



(5) 



where ^neutral Q i s defined by the expression 



t C/,D _ n/ U,Ds-i c U,DjrU,D 



1 In all next discussion we denote Neutral as £jv ,I? - 



(6) 
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Note that the charged couplings appear as a linear combinations of neutral couplings multiplied 
by Vckm matrix elements (more details see p3f). 

The next step is the calculation of the QCD corrections to b — > sl + l~ decay amplitude in the 
model III. The appropriate framework is that of an effective theory obtained by integrating out 
the heavy degrees of freedom. In the present case, t quark, W ± , H ± ) Hi, and H 2 bosons are the 
heavy degrees of freedom. Here H denote charged, Hi and H 2 denote neutral Higgs bosons. 
The QCD corrections are done through matching the full theory with the effective low energy 
theory at the high scale /i = mw and evaluating the Wilson coefficients from mw down to the 
lower scale /i ~ 0{m b ). Note that we choose the higher scale as /i = m w since the evaluation 
from the scale fi = mn± to fi = mw gives negligible contribution to the Wilson coefficients 
since the charged Higgs boson is heavy enough from the current theoretical restrictions, for 
example m H ± > 340 GeV |2§, m H ± > 480 GeV [f2§. 

The effective Hamiltonian relevant for our process is 



n 



eff 



-^v tb v t :Y,c^)o^) , (7) 



where the Oi are operators given in eqs. (||), @ and the Cj(/z) are Wilson coefficients renor- 
malized at the scale /i. The coefficients Ci(fi) are calculated perturbatively. 

The operator basis is the extension of the one used for the SM and 2HDM (model I and II) 
57] with the additional flipped chirality partners, similar to the ones used for the b — > S7 



decay in the model III [8J and SU(2) L x SU(2) R x 17(1) extensions of the SM @: 

01 = (s La 7 At C L( g)(c L( g7 M 6 La ), 

2 = {s L a r YfiCLa){cLp r y >J 'b L /3), 

q=u,d,s,c,b 
q=u,d,s,c,b 
q=u,d,s,c,b 



q=u,d,s,c,l 
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■fj,v 



°7 = -77r^s a a llu (m b R + m s L)b Q J : '' 

g _ 

09 = T^^La-ff.bL^ih^) , 
10TT Z 

g _ 



0x2 = (s La 'j^c La )(c R f 3 'y ll b R f 3 ), (8) 
and the second operator set 0[ — 0[ 2 which are flipped chirality partners of 0\ — 0\ 2 . 



0\ = 




0' 2 = 


{sRalfiCRa) {c R p^b R p), 


O's = 


q=u,d,s,c,b 


0' 4 = 


q=u,d,s,c,b 


oi = 


fsp^,^y,,6p^,) > ( a t R'y^' a t n) 

\ u Ha l/i'-'HaJ / < \HLp I HLp ) i 
q=u,d,s,c,b 


O'e = 


(sRal^Rp) (QL0Y<lLa}, 

n — 1 1 /i q f* hi 


0' 7 = 


\bix z 


O's = 


-f^SaT^a^mbL + m s R)b p g 


0' 9 = 


g _ 
1Q7T2 (sRal f ,bR. a )(l>y ll l) , 


o[ = 


g _ 
^^(sRalvbR^ih^l) , 


O'n = 


{sRonnC R p){cLf3^b L a) , 


0[ 2 = 


{sR n l i.CR a ){cLp r )' 1 b L p) , 



CL[IV 

1 



(9) 

where a and (3 are SU(3) colour indices and T^ v and Q^ u are the field strength tensors of the 
electromagnetic and strong interactions, respectively. 

Denoting the Wilson coefficients for the SM with Cf M (m\y) and the additional charged 
Higgs contribution with C^(mw), we have the initial values for the first set of operators (eq.(||)) 
ii 



C 'l,3,...6 > ll,12( m w) - 

C 2 3M (m w ) = 1 



c ' (mi " ) = K^W inx + 24(x-iy • 

C$ u (m w ) = — V«W + '"^"'' gW " D M + H' ■ 
szr^tV sin tV 9 

CgVwO = -r4-(5(x)-C(x)) , 



Cf... 6,11, 12 ( m w) 

C?(m w ) 
C»(m w ) 

Ci (m w ) 



+ 



+ 







mi 



■v. 



Is 



V 



tb 



V* 

r Ci 
V ts 



v. 



tb 



m t m h 

m t v ts v tb 

~ ~ + ^N,tcT^) (^N,bb + ^N,sbTr)^{y) > 
mt m b Vts V tb 



m 



V 



Is 



v tb 

~2 i^N,tt + £,N,tcT7^) (6iV",ti + ^N,tcTF~)I J \{y) 
"h v ts v tb 



and for the second set of operators eq. (|9]), 



Clf.6,ll,12( m H/) 

C' 7 H (m w ) 



C' s H (m w ) 
C' 9 H (m w ) 



. 
. 



2 (£w,fts777 + £,N,ss) (£,N,bb + ^sbTT - )^ 1 (^) 



+ 



Vb 
z D V tb 



'V 



tb 



(^AT,fes T A* + £.N,ss) (£iV,ti + £. 



AT,tc 

" tb 



V 



)F2{y) 



~2 KjV,6s TA * + £n,ss) (£.N,bb + £jV,s6 T/ )Gl(y) 
m t •'is •'tb 



+ 



V 

m t m b • '■'- J V ts 

z D V b 



(£N,bs^7; + €n,ss) (£n,U + £*N,tc 77^)^2(2/) 



m t v ts v tb 



D Vb 



m 



2 (£iV,6s 



+ £.N>s) (£,N,bb + ^N,sbT7~)Ll{y) 



ts 



'V 



tb 



where x = ml/m^y and y = mf/m 2 H ±. In eqs. (|10|) and (|TT1) we used the redefinition 



-vd 



I 4Gf w ,d 
V2 * 



The functions -B(x), C(x), D(x), i*i(2)(y), Gi(2)(y) ; H\(y) and are given as 



C(x) 
£>(ar) 



In x 



x — 1 (a: — 1)" 
x/2-3 + 3x/2 + 1 



lnx 



x — 1 (x — l) 2 
19x 3 /36 + 25x 2 /36 -x 4 /6 + 5x 3 /3 - 3x 2 + lGx/9 - 4/9 



(x-iy 



+ 



(x - l) z 



lnx 



Fi(y) 
F 2 (y) 
Gi(y) 
G 2 (y) 
Hi(y) 
Hy) 



y(7-5y-8y 2 ) y 2 (3y-2) 

72(2/ - l) 3 12(2/ -l) 4 Uy 
2/(52/ -3) , y(-3y + 2) 



12(y 



6(2/ -1)3 



In y , 



y(-y 2 + 5y + 2) 



+ 



24(2/ - I) 3 4(2/ 

y(y-3) y 

% - l) 2 + 2(2/ - 1) 
1 — 4sm 2 ^iy xy \ 1 
sin 2 9w 8 
1 £2/ I 1 



In 2/ 



In 2/ , 



1 



sin 2 9w 8 



2/-1 (2/-I) 2 
1 



In 2/ 



472/ 2 - 792/ + 38 32/ 3 - Qy + 4 



108(2/ -r 



18(2/ -i; 



In 2/ 



+ 



(2/-l) : 



In y 



(13) 



Note that we neglect the contributions due to the neutral Higgs bosons since their interactions 
include Yukawa couplings which should be very small (see the discussion section for details). 
For the initial values of the Wilson coefficients in the model III (eqs. (^)and (0)), we have 



2HDM ( \ 
l,3,...6,ll,m''WJ 


= o, 


Cl HDM {m w ) 


= 1 , 


C 2HDM (m w ) 


_ qSM 


Cl HDM {m w ) 


— U 8 


Cl HDM (m w ) 


r>SM 

— U 9 




— °10 



(m w ) + C 7 (m w ) 
(m w ) + Cg(m w ) 
(m w ) + Cg(m w ) 



■iH, 



I2HDM 
1,2,3,. ..6, 11, 12 


(m w ) 


= o, 








Q/2HDM 


(m w ) 


= C' 7 SM 


[m w ) 




[mw 


C1I2HDM 
°8 


(m w ) 


= c'£ M 


[m w ) 


+ c 8 H 


[mw 


r<l2HDM 
°9 


(m w ) 


= C' 9 SM 


[m w ) 


+ c 9 H 


[m w 


n l2HDM 
°10 


(m w ) 


riiSM 
— Wo 


[m w ) 


+ C[ H 


[m w 



(14) 

Using the initial values, we can calculate the coefficients C 2HDM and C' 2HDM at any lower 
scale with five quark effective theory where large logarithims can be summed using the renor- 
malization group and their evaluations are similar to the SM case. 

The operators 5 , 6 , On and 0\ 2 ( 0' 5 , 0' 6 , 0' u and 0' 12 ) give a contribution to the leading 
order matrix element of b — > 57 and the magnetic moment type coefficient C 7 (p) (C' 7 (/j,)) 



6 



is redefined in the NDR scheme as [28fl: 

C e 7 ff M = Cr DM {») + Q d (C 2HDM (») + N c C 2HDM (»)) , 

C 7 eff (») = C 7 2HDM (») + Q4C' 2HDM (») + N C C' 6 2HDM (»)) 

+ Qui— C? 2 HDM (») + N C ^ C' 2HDM (»)) . (15) 

TJlfj Tilly 

The NLO corrected coefficients C 2HDM (fi) and C' 2HDM (p) are given as 
C 2HDM ^) = C^ 2HDM ^) + ^-C? 2HDM ^) , 

c 7 2HDM {p) = c; LO ' 2 ™ f ( /U ) + ^c? 1)2HDM ( / i) . (i6) 

Arc 

The functions C 7 ,2HDM (n) and C' 7 L ,2HDM (/j 1 ) are the leading order QCD corrected Wilson 
coefficients [0, |2Bj, [Z7|: 

c Lo,2HDM {fi) = v ^c 2HDM (m w ) + (8/3)( V u / 23 - v ^)Ci HDM (m w ) 

i=i 

c ,LO,2HDM {fi) = ^6/23 C «HDM (mw) + (8/3) ^14/23_^6/23 )C «HDM (mw) 



where 77 = a s (mw)/o: s (fJ,), hi and cij are the numbers which appear during the evaluation [[Tql . 
C 7 ^ 2HDM (fi) is the a s correction to the leading order result that its explicit form can be found 
111 §|. C? 1)2 ^ A/ (/i) can be obtained by replacing the Wilson coefficients in Cj^ 2HDM (fi) 
with their primed counterparts. 

Since 2 (0 2 ) produce dilepton via virtual photon, its Wilson coefficient C 2 (/i) (C 2 (//)) 
and the coefficients Ci(/i), C^/z), C^/x) (C 2 (/i), C^(fi), C' 6 (fi)) induced by the operator 
mixing, give contributions to Cg^(fi) (Cg^(fi)). The perturbative part of Cg^(fi) fl5|, and 
Cg (n) including NLO QCD corrections are defined in the NDR scheme as: 

crv) = c 2 g HDM M(s) 

+ h(z, s) (3d(/i) + C 2 (jjl) + 3C 3 (/i) + C 4 (/i) + 3C 5 (/i) + C 6 (//)) 

- i/i(l,S)(4C 3 (/i)+4C74( / u) + 3C , 5(/i) + C , 6 ( / i)) (18) 

- h(0, s) (CM + 3C 4 (/i)) + 2 - (3C 3 (/x) + CM + 3C s (jm) + CM) , 

and 
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+ h( z , s) (3c;(/i) + cm + 3CM + cm + 3cm + cm) 

- ±h(lJ)(4CM+4CM + 3CM + CM) (19) 

- ~/i(o, s) (cm + 3C'M) + 1 (3C'M + cm + scM + CM) ■ 

2 

where z = — a-nd s = In the above expression 77(5) represents the one gluon correction 
to the matrix element O9 with m s = [[27] and the function h(z, s) arises from the one loop 
contributions of the four quark operators Oi,...,0 6 (0[, 0' 6 ). Their explicit forms can be 
found in fTTH WJ§ . In addition to the short distance part there exist the long distance (LD) effects 



due to the real cc in the intermediate states, i.e. the cascade process B — > K*ipi — > K*l + l~ 
where % = 1,..,6. In the current literature, there are four different approaches to take these 
intermediate states into account: HQET approach AMM approach |r2| |, LSW approach 
|32|1 and KS approach In our calculations, we follow the approach AMM. (We also use KS 



approach in our numerical analysis and see that results obtained in both methods are closed 
each other.) In this method the resonance cc contribution is parametrized using a Breit-Wigner 
shape with normalizations fixed by data, and this contribution is added to the perturbative one 
coming from the cc loop: 

C 9 e// (^) = C p 9 ert (fi) + Y reson (s) , (20) 
where Y reson (s) in NDR scheme is defined as 

Y (£) 3 V ^^"^ 



{3CM + c 2 (fj) + 3CM + cm + zcM + CM) ■ (21) 

For the expression C' 9 e ^(fi), it is enough to replace all unprimed coefficients with primed 
ones. In eqs. (18) and (19) the phenomenological parameter k = 2.3 is chosen to be able 
to reproduce the correct value of the branching ratio Br(B — > J/ipX — > Xll) = BriB — > 
J/ipX) Br(J/tp — > Xll) jnj. The NLO corrected coefficients Cj ,i — 1,...,6 can be found in 
p4| , |30]| . In our numerical calculations we neglect the coefficients C\ HDM (/x), C\ HDM (\i) and 
Cg 2 M (/i), Cq HDM (n) since they are numerically small at m& scale. 

Note that the expressions for unprimed Wilson coefficients in our case can be obtained from 
the results in |3(J by the following replacements: 

|y|2 _Lf7^ A_fU +^ —1 

l r I ^ 2 \^N,tt ' ^>N,tc ir t i \SJV,tt ^ ?iV,tc J ) 
m t ^ts K*6 
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m t m b V t * s ' Vtb 



To obtain primed coefficients, it is enough to use the primed ones at m\y level (eq. (|IID) 
since the evaluation of C-(fj,) from // = rriw to \i = ra^ is the same as that of Ci{p). 
For model II (model I) Y and XY are 

Y = l/tan/3 (l/tan/3) , 
XY = 1 (-l/tan 2 (3) . (23) 



At this stage we would like to make the following remark. In |M[ QED corrections to the 
Wilson coefficient for 5 — > X s ^ decay are calculated in addition to the NLO QCD ones 
and it was shown that QED corrections were almost one order smaller than NLO QCD ones. 
We expect that the similar situation exists for Cg . Therefore, in our calculations we neglect 
them. 

Finally, neglecting the strange quark mass, the matrix element for b — * s£ + £~ decay is 
obtained as: 

+ (C 10 (^)s7 m (1-7b)6 + ^ (/x)57 m (1+75)6) Hid (24) 

- 2fc 7 e// (//)^si<r^^ 
3 The exclusive B —> K*l + l~ decay 

Now, our aim is to look at the problem from the hadronic side. To calculate the decay width, 
branching ratio, etc., for the exclusive B — > K*l + l~ decay, we need the matrix elements 
(K* 1^7^(1 ± 7s)&| B), and (K* [sia^q^il ± 75)6! B), Using the parametrization of the form 
factors as in Jl9|], the matrix element of the B — > K*l + l~ decay is obtained as 



M = --^V tb V*\ hn \2A tot t pup(T ^p K ,q° + iB ltot e* - iB 2tot (e*q)(p B + p K ») p - iB 3tot (e*q)q^ 

2V Z7T [ L J 

+ Hid \zC tot e pupa e* u p K ,q a + %D xtot ^ - tD 2tot (e*q)(p B + p K *) p - iD ztot (e f q)q p ] j , (25) 

where e* M is the polarization vector of K* meson, p B and px* are four momentum vectors of B 
and K* mesons, q = Pb — Pk* and 

A tot = A + A' , 
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B\ tot 


= B 1 + B[ 


B2tot 


= B 2 + B' 2 


Bztot 


= B 3 + B' 3 


Ctot 


= c + c, 


Dim 


= D. + D', 


B>2tot 


= D 2 + D' 2 


Dstot 


= D 3 + D> 3 



-C e / f AC^T,, 

m B + m K * q 2 

-C e 9 ff (m B + m K .)Ai - 4C 7 e// ^(m| - m? K .)T 2 , 

-C^—^ AC^ ( T 2 + 2 ^ 2 T 3 

rriB + rriK* q \ m B — m K * , 



C e y-^f(A 3 -A ) + 4C 7 ^T 3 , 
V 



q 2 q 2 



10" 



-C 

mg + rriK* 
-C w (m B + m K *)A 1 , 

C M 
— wo ; , 

VtlB + THk* 

-C 10 ^(A 3 -A ) , 
q 2 



_ c >eff V 4 C^ / ^T 1 , 

m B + m K * q 2 

C' 9 eff (m B + m K *)A 1 + 4C; e// ^(m| - m 2 K .)T 2 , 

° 9 m im +4 ° 7 ^2"r 2 + ^2 Z^ 1 * 

m B -\- rriK* q \ m B — m K , / 

q 2 q 2 
V 



"ClQ 

ra s + rn K * 



C' 10 (m B + m K *)A 1 , 
Cl0 m B + m K , ' 

10 



D' 



C 



, 2m K * 



10 



(A 3 - Ao) 



(2f 



The hadronic formfactors V, Ai, A 2 , A , T 1; T 2 and T 3 has been calculated in the framework 



of light-cone QCD sum rules in j21], |22|.In our calculations we use the results of |22] where 
radiative corrections to the leading twist wave function and SU(3) breaking effects are also 
taken into account. The q 2 dependence of the form factors can be represented in terms of three 
parameters as: 



F(q 2 



ap- 



(29) 



where the values of parameters F(0),ap and bp are listed in Table 2. 





F(0) 


dp 


bp 


At 


0.34 ±0.05 


0.60 


-0.023 


A 2 


0.28 ±0.04 


1.18 


0.281 


V 


0.46 ±0.07 


1.55 


0.575 


Ti 


0.19 ±0.03 


1.59 


0.615 


T 2 


0.19 ±0.03 


0.49 


-0.241 


T 3 


0.13 ±0.02 


1.20 


0.098 



Table 1: The values of parameters existing in eq.(|29|) for the various form factors of the tran- 
sition B — > K*. 

Light cone QCD sum rules is applicable in the region, ml — q 2 ~ few GeV 2 and this 
leads to g^ aa; ~ 20 GeV 2 . To extend the results to the full physical region, we use the above 
parametrization in such a way that it correctly reproduces light cone QCD sum rules results 
up to q^ax ~ 20 GeV 2 . All the errors come from uncertainity of b-quark mass, higher twist 
wave functions, and variation of Borel paramater are added quadratically. 
Using eq.(p5|) we get the double differential decay rate: 
dT G 2 a 2 ^\V tb V t f\ 1/2 



dq 2 dz 



2 12 7r 5 m£ 



2Am| 



m 2 B s(l + z 2 )(\A tot \ 2 + \C tot \ 2 ) 



+ 



+ 



2r 
1 

2r 



Am|(l -z 2 ) 



B 2t ot\ + \D 2 tot\ 



'l tot 



- z 2 ) + 8rs} (\B ltot \ 2 + \D, 
2 \ m %(l -r-s)(l- z 2 ) {Re (B ltot B* 2tot ) + Re (D ltot D* 2tot )} 
{Re(B ltot C; ot )+Re(A tot D* ltot )} 



(30) 
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where z = cosO , 6 is the angle between the momentum of £ lepton and that of B meson in the 



center of mass frame of the lepton pair, A = 1 + r 2 + s 2 — 2r — 2s — 2rs, r = and s 



m, 



Q 

m. 



As we noted, Afb and Pi can give more precise information about the Wilson coefficients 
C 7 e// , C 9 e// and Cio. They are defined as: 



i dT 

a dq 2 dz 



dT 

-i dq 2 dz 



1 dz dT 
ii dq 2 dz 



-i dq 2 dz 



(31) 



and 



-i: 



fie/ 



2« 



(32) 



■1^ 



rir _ ^ fir 

fig 2 fig 2 

where £ = — 1(+1) corresponds to the left (right) handed lepton in the final state. After the 
standard calculation, we get 



1 f 32 

A jy 

4 Xm% 



s\m%Re (A tot C; ot ) + \^-Re (R 



2totD* 2tot )+A-^ 



A 



+ Ars 



1 _ r _ s ) [jRe (5 2totJ D* t J + i?e (S ltot £> 2 * tot )] 



ife (S UotJ D* tot ) 
(33) 



The denominator A in eq. (^) can be obtained from eq.(|3"UD by integration over z of the 
terms within the curly bracket. In our numerical analysis we used the input values given in 
Table ©. 



Parameter 


Value 


m c 


1.4 (GeV) 


m b 


4.8 (GeV) 


a- 1 

"em 


129 


At 


0.04 




3.96 • 10~ 13 (GeV) 


m Bd 


5.28 (GeV) 


m t 


175 (GeV) 


m w 


80.26 (GeV) 


m z 


91.19 (GeV) 


h-QCD 


0.214 (GeV) 


a s (m z ) 


0.117 


sinOw 


0.2325 



Table 2: The values of the input parameters used in the numerical calculations. 
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4 Discussion 



Before we present our numerical results, we would like to discuss briefly the free parameters of 
the model III. This model induces many free parameters, such as t^- D where i,j are flavor indices 
and it is necessary to restrict them using the experimental measurements. The contributions 
of the neutral Higgs bosons Hq and A to the Wilson coefficient Cj (see the appendix of [3~7] for 
details) are 



D 



l= d,s,b ' 8 numb 



Cf°{m w ) = {V tb V t :r l £ % % -9*—, (34) 

^ b 8m imb 

where rrti and Qi are the masses and charges of the down quarks (i = d, s, b) respectively. 
These expressions show that the the neutral Higgs bosons can give a large contribution to the 
coefficient C-j which is in contradiction with the CLEO data announced recently fl38|, 



Br(B -> X s7 ) = (3.15 ± 0.35 ± 0.32) 10" 4 . (35) 

Such dangerous terms are removed with the assumption that the couplings ^ is {i = d, s, b) and 
£jvd& are negligible to be able to reach the conditions ^N,bb^N,is << 1 anc i ^N,db^N,ds << 1- With 
the discussion given above and using the constraints [pE|j , coming from the AF = 2 mixing, the 
p parameter and the measurement by CLEO Collaboration we get following restrictions: 
£mc << ^NttJ ^Nbb an d £mb ~ , - ~ 0, where the indices i, j denote d and s quarks . Under 
this assumption, the Wilson coefficients C' 7 , C' g and C' 1Q can be neglected compared to unprimed 
ones and the neutral Higgs contributions are supressed because the Yukawa vertices are the 
combinations of £jy s6 and £jv ss . After these preliminary remarks, let us start with our numerical 
analysis. 

In this section, we study the q 2 dependencies of the differential Br, A FB and Pi of the 
decay B — > K*l + l~ for the selected parameters of the model III Cmb)- I n fig s - an d 

|2] we plot the differential Br of the decay B — > K*l + l~ with respect to the dilepton mass 
q 2 for £x bb = 40 m b and charged Higgs mass m H ± = 400 GeV at the scale /x = m b . Fig. |l] 
represents the case where the ratio \r t b\ = \jff^\ « 1- I n this case, it is shown that the 
differential Br obtained in the model III almost coincides with the one calculated in the SM. 
In Fig. 0, we present the same dependence as in Fig. [1| for the case r tb » 1 and ^ bb = 40 m b . 
Here ^% tt lies in the restriction region coming from the constraints mentioned above. The 
strong enhancement is observed compared to the SM result especially for the small values of 
q 2 . Further, the differential Br increases with the increasing ^Nbb f° r r tb >> 1- 
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Figs. ^| and |] show the q 2 dependence of A FB for ^ bb = 40 m b and charged Higgs mass 
m#± = 400 GeV at the scale /i = m b . For \r t b\ « 1 (Fig. |3|) Afb changes its sign, however 
for rtb >> 1 (Fig. (|) it is positive withouth LD effects. Therefore the determination of the 
sign of Aps in the region < q 2 < 2.8 (GeV) 2 (here the upper limit corresponds to the value 
where Afb change sign in the SM) can give a unique information about the existence of the 
model IIL Further, Afb becomes more positive with increasing ^ bb . It is interesting to note 
that the final lepton tends to move almost in the direction of B meson with increasing ^ bb for 
r t b » 1, especially in the q 2 region far from LD effects. 

The behaviour of Pl versus q 2 is presented in series of figures (|| and H) f° r £,N,bb an d charged 
Higgs mass = 400 GeV at the scale fi = m b . Pl becomes more negative in the region 

r tb >> 1 (Fig |B|) compared the one in the region \r tb \ << 1 (Fig ^). This tendency increases 
with increasing £,® bb . 

Now we would like to present the values of Br for the B — > K*l + l~ decay in the SM and 
model III, without the LD effects. After integrating over q 2 , we get 

Br(B -> K*l + l~) = 0.278 ± 0.050 x 10~ 5 (SM) (36) 

and for the model III 

f 0.278 ± 0.050 x 10~ 5 (\r tb \ « 1) 
Br(B -> K*l + l~) = I (37) 

[ 0.990 ± 0.050 x 10~ 5 (^ » 1, ^ 6b = 40m 6 ). 

The Br in the SM is almost equal to the model III value for \r tb \ « 1. However, the strong 
enhancement can be observed for r t b » 1, especially with increasing ^ bb . Using the CLEO 
data for the inclusive branching ratio Br(B — > X s l + l~) < 0.42 10~ 4 [H(J, it is possible to 
estimate an upper limit for £j} bb , which is £,®,bb — 40 m b . 

In conclusion, we analyse the selected model III parameters ( ^§ bb , £%f U ) dependencies 
of the differential Br ,Ap B and Pl of the decay B — > K*l + l~. We obtain that the strong 
enhancement of the differential Br is possible in the framework of the model III and observe 
that Aps and Pl are very sensitive to the model III parameters (^ 66 , Cnu)- Therefore, their 
experimental investigations ensure a crucial test for new physics. 
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10*7 dBr/dq A 2 





q A 2 (GeV A 2) 



Figure 1: Differential Br as a function of q 2 for fixed £§ bb = 40 m b in the region \r tb \ << 1, at 
the scale fx = m b . Here solid line corresponds to the model III with LD effects, and dashed line 
to the SM withouth LD effects. 



10*7 clBr/dq A 2 
10 





q"2 (GeV A 2) 



Figure 2: The same as Fig 1, but at the region r tb » 1. Dotted dashed line corresponds to 
the model III withouth LD effects. 
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Figure 6: Pl as a function of q 2 for fixed £^ b6 = 40 in the region r t b » 1, at the scale 
11 = 171},. Here solid line corresponds to the model III with LD effects, and dashed line to the 
model III withouth LD effects. 
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